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PORTFOLIO STRUCTURING USING 
LOW -DISCREPANCY DETERMINISTIC SEQUENCES 

The United States Government has certain rights to this 
invention pursuant to grants CCR-91-14 042 and IRI- 92 -12597 
awarded by the National Science Foundation, and to grant 
AFOSR-91-0347 awarded by the U.S. Air Force. 
5 This application is a continuation-in-part of pending 

0^ application No. 08/902 , 921 J/v which was filed July 30, 1997 as 

a continuation application of now-abandoned application No. 
08/285,902 filed August 4, 1994. 



fti 10 Background of the Invention 

£j The invention relates to portfolio structuring and, 

: „ i 

r{ more particularly, to structuring of a portfolio including 

* assets/liabilities such as real estate, stocks, bonds, 

jij futures, options, collateralized mortgage obligations, 

*% 15 single-premium deferred annuities and insurance contracts. 

Portfolio structuring includes transactions such as 
buying, holding and selling of assets/liabilities, as well 
the adoption, continuation, modification or termination of 
transaction strategies. Typically, decisions concerning 
2 0 transactions and strategies are based on the evaluation of 
one or several complicated functions of a large number of 
variables such as a mult i -dimensional integral or the 
inverse of a distribution function, for example. 

In financial securities transactions, which includes 
25 the initial sale, the value of a security may be estimated, 
e.g., based on expected future cash flow. Such cash flow 
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may depend on variable interest rates, for example, and 
these and other relevant variables may be viewed as 
stochastic variables. 

It is well known that the value of a financial security 
which depends on stochastic variables can be estimated in 
terms of a mult i -dimensional integral. , The dimension of 
such an integral may be very high. 

In collateralized mortgage obligations (CMO) , for 
example, instruments or securities variously called 
tranches, shares, participations, classes or contracts have 
cash flows which are determined by dividing and distributing 
the cash flow of an underlying collection or pool of 
mortgages on a monthly basis according to pre- specif ied 
rules. The present value of a tranche can be estimated on 
the basis of the expected monthly cash flows over the 
remaining term of the tranche, and an estimate of the 
present value of a tranche can be represented as a multi- 
dimensional integral whose dimension is the number of 
payment periods of the tranche. For a typical instrument 
with a 3 0 -year term and with monthly payments, this 
dimension is 360. 

Usually, such a high-dimensional integral can be 
evaluated only approximately, by numerical integration. 
This involves the generation of points in the domain of 
integration, evaluating or "sampling" the integrand at the 
generated points, and combining the resulting integrand 
values, e.g., by averaging. Well known for numerical 
integration in securities valuation is the so-called Monte 
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Carlo method in which points in the domain of integration 
are generated at random. 

With integrands arising in financial securities 
valuation, the computational work in combining the sampled 
values is negligible as compared with producing the 
integrand values. Thus, numerical integration methods in 
securities valuation may be compared based on the number of 
samples required for obtaining a sufficiently accurate 
approximation to the integral. 

Sampling techniques are useful also as applied to 
formulations which need not be in terms of an integral. In 
particular, in portfolio structuring, this applies to a 
quantity known as value at risk (VAR) which may be defined 
as the worst potential loss of value of a portfolio over a 
period of time under normal market conditions, at a 
specified confidence level. While value at risk can be 
expressed as a definite limit in a one-sided indefinite 
integral (see Philippe Jorion, Value at Risk , Irwin 
Professional Publishing, 1997, p. 88) , a more practical 
formulation is in terms of the inverse of a cumulative 
distribution function. 

Summary Of The Invention 

A preferred method for estimating the value of 
a financial security involves numerical integration unlike 
Monte Carlo integration in that an integrand is sampled at 
deterministic points having a low-discrepancy property. As 
compared with the Monte Carlo method, significant advantages 
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are realized with respect to speed, accuracy, and 
dependabi 1 i t y . 

A computed estimate can be used in setting an offering 
price, e.g., somewhat higher that the estimate. An estimate 
5 can be used also in deciding whether to buy, hold or sell 
the security. For example, if the estimate exceeds the 
price at which the security is being offered for sale, a 
"buy" recommendation is appropriate. 

A computed estimate can represent the value of other 
10 complex securities, e.g. a debt or insurance instrument. 

Further to estimating the value of a complex security, 
sampling at points of a low-discrepancy deterministic 
sequence can be used in estimating value at risk in 
portfolio structuring . 

15 

Brief Description Of The Drawing 

Fig. 1 is a schematic of a programmed computer system 
in accordance with a preferred embodiment of the invention. 

Fig. 2 is a graphic representation of performance data 
2 0 obtained in computer trial runs with an exemplary embodiment 
of the invention as compared with two Monte Carlo 
computations . 

Fig. 3 is a flow diagram of a preferred embodiment of a 
method in accordance with the invention. 
25 Fig. 4 is a flow diagram of a preferred further 

embodiment of the invention, including the use of multiple 
processors. 
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Fig. 5 is a flow diagram of initially selling a complex 
security in accordance with a preferred embodiment of the 
invention. 

Fig. 6 is a flow diagram of buying/holding/selling a 
5 complex security in accordance with a preferred further 
embodiment of the invention. 

Fig. 7 is a flow diagram of determining value at risk 
(VAR) for portfolio structuring in accordance with a 
preferred embodiment of the invention. 
10 Fig. 8 is a flow diagram of an example of portfolio 

structuring in accordance with a preferred embodiment of the 
invention. 

Further included is an Appendix with two computer 
algorithms in "C" source language, respectively for 
15 computing Sobol points and Halton points. For a description 
of C, see B.W. Kernighan et al . , The Programming Language C , 
Prentice -Hall, 1978 . 

Detailed Description Of Preferred Embodiments 
2 0 Fig. 1 shows a stored-program computer 11 connected to 

input means 12, e.g., a keyboard, for entering financial 
securities data, and connected to output means 13, e.g., a 
visual display device, for displaying an estimated value of 
the financial security. The computer 11 includes a working 
25 memory 111 (M) , a low-discrepancy deterministic point 

generator 112 (P) , an integrand evaluator 113 (E) , and an 
integrand-value combiner 114 (C) . 

In estimating the value of a mult i -dimensional integral 
in financial securities valuation, a multivariate integrand 
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is sampled at points corresponding to a low-discrepancy 
deterministic sequence of points in the multivariate unit 
cube as defined below. If the multivariate unit cube is 
also the domain of integration, the points of the low- 
discrepancy deterministic sequence can be used as sample 
points directly. In the case of a more general domain of 
integration, sample points correspond to points of a low- 
discrepancy deterministic sequence in the multivariate unit 
cube via a suitable transformation or mapping. 

When a sufficiently large number of sampled values has 
been computed, an approximation of the integral is obtained 
by suitably combining the computed values, e.g., by 
averaging or weighted averaging. 

In the d-dimensional unit cube D = [0,l] d , a low- 
discrepancy deterministic sequence z 17 z 2 , ... of points 
in D can be characterized as follows: 

For a point t = [t lt . . . , t d ] in D, define 
[0, t) = [0, t 2 ) x . . . x [0, t d ) , 
where [0, t ± ) denotes a left-closed, right-open interval, 
and denote with Xto, t) ( • ) the characteristic or indicator 
function of [0, t) . For points z lf . . . , z n in D, define 

Rn (t /" Z x , . . . , Z n ) = (£koi n X[0, t) ( Z k) ) / n ~ t x t 2 . . - t d , 

and define the discrepancy of z lt . . . , z n as the L ro -norm of 
the function R n ( . ; z lf . .., z n ) , i.e., 

||R n (.; z x , z n )|| TO = sup teD |R n (t; z lf z n ) | . 

The sequence z lt z 2 , ... is said to be a low- 
discrepancy deterministic sequence provided 

||R n (. ; z lf z n ) |U - 0( (log n) d /n) . 
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Low- discrepancy deterministic sequences are described 
in the published literature; see, e.g., H. Niederreiter , 
"Random Number Generation and Quasi-Monte Carlo Methods", 
CBMS-NSF , 63, SIAM, Philadelphia, 1992 and S. Paskov, 
5 "Average Case Complexity of Multivariate Integration for 

Smooth Functions", Journal of Complexity , Vol. 9 (1993), pp. 
291-312. Well-known examples of low-discrepancy 
deterministic sequences are the so-called Hammersley points, 
Halton points, Sobol points, and hyperbolic-cross points. 
10 For illustration, in the case of Sobol points in a 

single dimension (d=l) , a constructive definition may be 
given as follows: Choose a primitive polynomial 

P(x) = x n + aiX" -1 + . . . + a n .iX + 1 
(whose coefficients a ± are either 0 or 1) and define so- 
15 called direction numbers v if i > n by the following 
recurrence formula : 

v ± = a{v ± . x © a 2 Vi_ 2 © . . . © a^v^! © Vi_ n © (Vi_ n /2 n ) , 
where © denotes a bit -by-bit "exclusive or" operation. 

The initial numbers v x = n^/2, . . . , v n = m n /2 n can be 
2 0 chosen freely provided m ± is odd and 0 < mi < 2 1 for i = 
1 , 2 , . * . , in . 

Using the direction numbers v ± so defined, now define 
the one-dimensional Sobol sequence x x/ x 2 , ... by 
x k = b x v x © b 2 v 2 © ... © b w v w/ k ^ 0 
25 where k = £i=o ri ° 9 k1 b i 2 i is the binary representation of k. 

For higher dimensions (d > 1) , the first d primitive 
polynomials P lf P 2 , P d are used. If {xj/}^ 00 denotes 

the one-dimensional Sobol sequence generated by the 
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polynomial P i; the d-dimensional Sobol points are defined as 

x k = ( X k X / X k 2 / • • • / x k d ) • 

While this definition can be implemented as a computer 
algorithm directly, faster techniques are known which 
produce these points in a "shuffled" or permuted sequence. 
In particular, this applies to the computer algorithm given 
in the Appendix. 

For specificity in the following, a preferred 
embodiment of the invention is described as applied to a 
collateralized mortgage obligation known as CMO FN, 89-23. 
This has thirty-year maturity and consists of the following 
tranches : 

PAC tranches 23-A, 23-B, 23-C, 23-D, 23-E 

supporting tranches 23 -G, 23 -H, 23 -J 

residual tranche 23 -R 

accrual tranche 23 -Z 
The monthly cash flow is divided and distributed according 
to pre-specif ied rules which are included in a formal 
prospectus. Some of the basic rules may be stated as 
follows : 

First from the monthly cash flow, the coupon is paid to 
the tranches. The remaining amount, called Principal 
Distribution Amount, is used for repayment of the principal. 
Prior to a fixed date in the future, the Principal Amount 
will be allocated sequentially to the tranches 23 -G, 23 -H, 
23- J and 23-Z. After that date, the Principal Distribution 
Amount will be allocated sequentially to the tranches 23-A, 
23-B, 23-C, 23-D and 23-E according to a planned schedule. 
Any excess amount of the Principal Distribution Amount over 
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the planned schedule will be allocated sequentially to the 
tranches 23-G, 23-H, 23- J and 23-Z. A distribution of 
principal of the tranche 23 -R will be made only after all 
other tranches have been retired. 
5 In deriving an estimate for the present value of the 

security at the time of issue, the following notation is 
used below: 

C - the monthly payment on the underlying mortgage 

pool ; 

10 i k - the projected interest rate in month k, k = 1, 2, 

P . . . , 360; 

:rj w k the percentage of mortgages prepaying in month k; 

^ a 360 - k+ i " t ^ ie remaining annuity after month k. 

U 

%J A remaining annuity a k can be expressed as 

y 15 a k = 1 + v 0 + . . . + v 0 k_1 , k = l, 2, . .., 360, 

JL with v 0 = l/(l + i 0 ), where i 0 is the current monthly 

RJ interest rate. Thus, after k months, the remaining amount 

of principal borrowed is C-a k . 
3S" It is assumed that the interest rate i k can be 

20 expressed as 

i k = K 0 exp(^) i t _ lt 
where exp ( . ) denotes exponentiation and where 
£,!, ^ 2 , ^ 360 are independent, normally distributed random 

variables with mean 0 and variance a 2 , and K 0 is a given 
25 constant. For the present example, a 2 = 0.0004 is chosen. 

It is assumed further that w k as a function of i k can 
be computed as 

w k = K x + K 2 arctan (K 3 i k + K 4 ) , 
where K lf K 2 , K 3 , K 4 are given constants. 
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Under these assumptions, the cash flow in month k, k = 
1 , 2 , . . . , 360 is 

Cd-WitU) • • ■ (l-w^^i, • - - ,5 k -i)) (l-w k (? 1# . - ., W 

5 where 

w k Ui* = Ki + K 2 arctan(K 3 K 0 k i 0 exp(^ 1 +. . .+£ k ) +K 4 ) . 

This cash flow is distributed according to the rules of 
FN, 89-23. Then, the cash flow for each of the tranches is 
multiplied by the discount factor 
10 v x . . -v k (^, . . . , £ k ) , with 

Vj . . . , Sj) = l/d + K 0 U 0 exp(5 1 +. . .+£.,) ) , j=l,2, . . . ,360, 
to find the present value for month k. Summing of the 
present values for every month gives the present value PV X , 
for each tranche T. 
15 The expected or estimated value, 

E(PV T ) = E(PV T (^, . . .,5 36 o>> , 
upon a change of variables is represented by 

E(PV T ) = J D PVTfy^Xi) , . . . ,y 36 o(x3 60 ) )dx x . . .dx 360 / 

where 

20 Xi = (200)"^./ exp(-t 2 / (2a) ) dt . 

Thus, for each tranche T, a 360-variate integrand has 
to be integrated over the 3 60 -dimensional unit cube, . 
After generating a point 

(x-l , x 2 , . . . , x 360 ) 
2 5 of a low-discrepancy deterministic sequence in the unit 
cube, the point 

(Yi/ Y2, • • - , Y360) 
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is computed by finding the value of the inverse normal 
cumulative distribution function at each x k . Then, for each 
tranche T, the function value 

pv t (yi, y 2 / • • • / Y360) 

is computed. These are the function values used in 
numerical integration . 

Fig. 2 shows results from trial runs for CMO FN, 89-23 
with a preferred method using Sobol points generated by the 
corresponding computer algorithm given in the Appendix, as 
compared with Monte Carlo integration. Two Monte Carlo 
computations were carried out, with different "seeds" or 
starting values of a congruential pseudo-random number 
generator known as RAN2 ; see W. Press et al . , Numerical 
Recipes in C , Cambridge University Press, 1992. It is 
apparent that the preferred method reaches a steady value 
more rapidly. In this and other trial runs, with typical 
requirements of precision and confidence, a speed-up by a 
factor of 3 to 5 was realized as compared with Monte Carlo 
integration. Much greater speed-up can be expected when 
higher precision or/and higher levels of confidence are 
sought . 

In a further trial run with CMO FN, 89-23, instead of 
Sobol points, Halton points were used as generated by the 
corresponding computer algorithm given in the Appendix. It 
is felt that Sobol points may be preferred over Halton 
points for integrals of high dimension. However, this 
preference may not apply in the case of lower-dimensional 
integrals, e.g., with dimension up to 5 or so. 
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A computation as described may be terminated after a 
predetermined number of function evaluations. 
Alternatively, e.g., after every function evaluation or 
after a predetermined incremental number of function 
evaluations, a current approximation may be compared with 
one or several preceding approximations, for termination 
once a suitable condition depending on the difference 
between approximations is met. Such termination criteria 
may be called automatic. Automatic termination is 
particularly reliable where a sequence of approximations 
settles down smoothly; see, e.g., the curve in Fig. 2 
corresponding to Sobol points. 

Advantageously in computing function values, a cluster 
or network of multiple parallel processors or workstations 
can be used. This may involve a master or host processor 
providing points of a low-discrepancy sequence to slave 
processors and combining function values returned by the 
slave processors into an approximate value for the integral. 
Thus, the computation can be speeded up in proportion to the 
number of processors used. 

Advantageous further, in combination with numerical 
integration as described above, is the use of error 
reduction techniques analogous to variance reduction in 
Monte Carlo integration as described, e.g., by M. Kalos et 
al . , Monte Carlo Methods , John Wiley & Sons, 1986. This may 
involve a change of variables or/and variation reduction, 
for example. 

The specific embodiment described above, involving 
valuation of a collateralized mortgage obligation, can be 
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interpreted as follows: The definite integral represents 
the expected value of the sum of a series of payments on the 
obligation, namely the regular payments of mortgage 
principal and interest, plus any mortgage prepayments. The 
likelihood of prepayment depends in part on the prevailing 
level of interest for mortgage refunding. 

For each payment, the value of the payment is weighted 
by the probability of all interest rate scenarios up to the 
time of the payment. Thus, for example, the amount of next 
month's payment will be valued based on the probability of 
all possible interest rates likely to prevail next month, 
and the impact of each of those possibilities on the 
likelihood of prepayment of principal. For the month after, 
valuation takes into account the interest rate scenarios of 
both months, their probabilities of occurrence, and the 
effect of interest rates on prepayment. In addition to next 
month's valuation, the month-after estimate will include an 
additional term corresponding to the month-after payment, 
represented by an integral over the possibilities of the 
two-month process. Further monthly payments, i.e. third, 
fourth et seq., are represented analogously. With each 
additional integral term representing a scheduled or 
unscheduled payment amount, the over-all value of the CMO is 
represented by the sum of the terms . 

The value of the CMO can be decomposed further, e.g., 
by valuing the interest payment and the principal payment 
separately in each period to represent an "interest -only 
tranche" and a "principal -only tranche". Thus, the value of 
the CMO can be represented by the sum of a yet -greater 
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number of integrals. Conveniently, the value of each of the 
individual integrals in the sum can be estimated as 
described above, using a low-discrepancy deterministic 
sequence of sample points. 

A sum of payment integrals can also be used to 
represent the value of a portfolio of two or more CMOs, 
giving the same result as a single integral over the sum of 
the payments. The range of possibilities over which the 
integration is effected is the same in both cases. 

A CMO portfolio can be augmented in other ways, e.g., 
by the inclusion of one or several bonds. For valuation 
purposes, a high-grade bond differs from a CMO in that 
prepayment will take place for the entire issue or not at 
all. Again, valuation depends on the fluctuation of 
interest rates, taking into account all possible future 
interest rate scenarios. This applies whether there is one 
payment (from a single CMO) , several payments (from several 
CMOs) , or payments from a combination of CMOs and bonds. 

Other than CMOs and bonds, a complex security portfolio 
may include securities such as single-premium deferred 
annuities (SPDA) . An SPDA is a contract sold by an 
insurance company to an annuitant investor. Under the 
contract, the insurance company takes money from the 
investor, invests this money in a portfolio of various 
investments, and credits a portion of the investment 
proceeds to the investor's account. The SPDA is an asset of 
the insurance company. There is an initial investment by 
the company in the cost of producing the contract. There is 
a predictable return for the company on the investment 
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portfolio, of interest and principal payments by the various 
borrowers less the contractual crediting to the annuitant's 
account. The unpredictable factor in the valuation of the 
portfolio is again random movements of interest rates, as 
5 these influence prepayment of investment assets as well as 
unexpected withdrawals and cash demands by the annuitant. 
The same methodology of valuation can be used after 
inclusion of the SPDA asset in the complex security. 

Valuation, by the insurance company, of a portfolio 

10 including an SPDA contract involves estimating the value of 
the contract itself. In pricing a new contract, its value 
is estimated in advance of issue by making assumptions about 
the kinds of assets that will be in the portfolio along with 
the SPDA. The pricing of the SPDA before issue and its 

15 valuation after issue are equivalent to the calculation 
presented above for a CMO. 

A life insurance policy, or any other kind of insurance 
contract differs from the above in that, in addition to the 
interest rate variable, there is an additional stochastic 

20 variable representing health risk, mortality risk or fire 
risk, for example. Under inclusion of this variable, the 
same techniques as described above can be used for 
valuation. 

Many other types of securities are amenable to 
25 corresponding valuation and estimation, where each period's 
prospective payment is valued depending upon all previous 
values of a stochastic variable. For example, in the case 
of an interest rate option, the stochastic variable is 
interest rate and the payment is the value of a bond. In 
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the case of a stock option, the variable is the price of the 
stock. In the case of a currency option, the variable is 
the relative value of two currencies. 

As an example of an application of sampling of a 
5 function other than an integrand, the following concerns 
value at risk (VAR) . If 

h = h(5i, . . . , U 
denotes the potential loss of a portfolio whose value 
depends on d market factors represented by stochastic 

10 variables £ 1# . .., £ d/ and if F h denotes the cumulative 
distribution function of h, i.e. 

F h (x) = Prob {h < x} , 
then, for a specified percentile or quantile p, 

VAR = x p = F h _1 (p) . 

15 The value of p represents a level of confidence, with p=l 

representing certainty. For practical purposes, a choice of 
p=0.98 is typical, indicating that the loss may exceed VAR 
with a probability of 2%. 

In accordance with a preferred embodiment of the 

20 invention, VAR can be estimated by a technique as 

illustrated by Fig. 7, using sampling of the potential loss 
function at m points of a low-discrepancy deterministic 
sequence. The potential loss h k is evaluated for each point 
of the sequence, resulting in a list h lf h m . The 

25 number of points, can be pre-set, e.g. at m=1000. The list 
is sorted, and the cumulative distribution function is 
determined. The cumulative distribution fucntion is a 
monotonically increasing function whose domain {x} includes 
the h-values and whose range {F h (x) } is in the interval from 
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0 to 1 . The desired estimate of VAR is determined as that 
value of x for which F h (x)=p. (It is understood that, 
wherever exact equality cannot be realized in calculation, a 
value will be approximated or interpolated.) 
5 An implementation of the technique was applied to a 34- 

dimensional problem, namely a portfolio of 34 at-the-money 
currency and equity options having maturity of one year, 
with $10,000 allocated to each of the options, yielding a 
total portfolio of $340,000. Volatility and correlation 
10 data were taken from J.P.Morgan/Reuters RiskMetrics data 
sets . 

The technique was applied further to a 360 -dimensional 
problem, namely the CMO portfolio consisting of ten tranches 
as described above, with $100,000 allocated to each tranche, 

15 yielding a total portfolio or $1,000,000. 

For the options portfolio, 1,000 points of the low- 
discrepancy deterministic generalized Faure sequence 
resulted in VAR being determined to within 1% of a reference 
determination using 100,000 points. For the CMO portfolio, 

20 a determination of VAR using 1,500 points of the same 

sequence also was within 1% of a reference value, obtained 
from 150,000 points. These results are contrasted with 
prior-art Monte-Carlo computations where, for the same 
numbers of samples (1,000 and 1,500, respectively), VAR 

25 could be estimated only to within 10% of the respective 

reference values, i.e. with 10 -fold error as compared with 
the low-discrepancy deterministic sequence determinations. 

In the technique as described, computational effort is 
greatest in evaluating the h k 's, and the computation can be 
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accelerated significantly by using multiple processors for 
evaluation in a fashion analogous to Fig. 4. In this 
fashion, e.g. if 10 processors are available, a 10-fold 
increase can be achieved in the computational speed in 
sampling of the loss as compared with a single processor. 

As illustrated by Fig. 8, a computed value of VAR can 
be compared with a target value which is considered 
acceptable. If VAR meets the target level, the portfolio 
may remain unchanged; if VAR is considerably less than the 
target level, the portfolio may be restructured for 
assumption of greater risk and potentially increased return; 
if VAR exceeds the target level, restructuring for reduced 
risk is indicated. Restructuring may be performed 
immediately or over time, e.g. by changing investment 
policies affecting future transactions. 
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